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ABSTRACT 

We select 25,000 galaxies from the NEWFIRM Medium Band Survey (NMBS) to study the rest- 
frame U — V color distribution of galaxies at < z < 2.5. The five unique NIR filters of the NMBS 
enable the precise measurement of photometric redshifts and rest-frame colors for 9,900 galaxies at 
1 < z < 2.5. The rest-frame U — V color distribution at all z < 2.5 is bimodal, with a red peak, a 
blue peak, and a population of galaxies in between (the green valley). Model fits to the optical-NIR 
SEDs and the distribution of MlPS-detected galaxies indicate that the colors of galaxies in the green 
valley are determined largely by the amount of reddening by dust. This result does not support the 
simplest interpretation of green valley objects as a transition from blue star-forming to red quiescent 
galaxies. We show that correcting the rest-frame colors for dust reddening allows a remarkably clean 
separation between the red and blue sequences up to z ^ 2.5. Our study confirms that dusty starburst 
galaxies can contribute a significant fraction to red sequence samples selected on the basis of a single 
rest- frame color (i.e. U — V), so extra care must be taken if samples of truly "red and dead" 
galaxies are desired. Interestingly, of galaxies detected at 24 /im, 14% remain on the red sequence 
after applying the reddening correction. 

Subject headings: galaxies: formation — galaxies: evolution — galaxies: high-redshift 



1. INTRODUCTION 

A powerful and observationally inexpensive tool for 
studying galaxy evolution is the color-magnitude or, 
more fundamentally, the color-mass diagram, in which 
the distribution of galaxies is bimodal with early-types 
residing on a well-defined sequence at red colors that 
is largely separate from a cloud of blue late-types (e.g. 
Blanton et al. 2003). Recent deep surveys have enabled 
the study of the evolution over < z < 1 of luminos- 
ity (Bell et al., 2004, hereafter B04; Brown et al. 2007; 
Faber et al. 2007) and mass (Borch et al. 2006) functions 
divided between red and blue samples selected based on 
color-magnitude and color-mass diagrams. A primary re- 
sult of these studies is that the stellar mass in red galaxies 
builds up by a factor of ~2 between z = 1 and z = 0, 
while the stellar mass in blue galaxies remains roughly 
constant over the same redshift range. This may imply 
a migration from the blue to red sequences (B04; Faber 
et al. 2007). In the simplest models the evolution from 
blue to red colors is largely driven by simple aging of stel- 
lar populations, with the rate of objects entering the red 
sequence modulated by the particular feedback processes 
at work that regulate star formation (e.g., Schawinski 
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et al. 2007). 

Recently, several surveys have started to push the sys- 
tematic study of galaxy colors to z > 1. There is 
some evidence that the color bimodality persists up to 
z = 2, though the contrast separating the red and blue 
sequences is low. This is either because spectroscopic 
samples that span a broad range of colors are small (Gi- 
allongo et al. 2005; Franzetti et al. 2007; Cassata et al. 
2008; Kriek et al. 2008) or because larger photomet- 
ric samples suffer from uncertain photometric redshifts 
at z > 1.5, which leads to uncertain rest-frame colors 
(Taylor et al. 2009a). Wuyts et al. (2007) (hereafter 
W07) and WiUiams et al. (2009) (hereafter W09) demon- 
strate that quiescent galaxies up to z ~ 2 can be fairly 
cleanly separated from the actively star-forming popula- 
tion when multiple rest-frame colors are used. 

In this Letter we use the NEWFIRM Medium Band 
Survey (NMBS; van Dokkum et al. 2009) to explore the 
evolution of the galaxy color bimodality over < z < 2.5. 
Wc show that two independent methods for accounting 
for dust in objects with intermediate colors allows the 
separation of a clearly-defined red sequence up to at least 
z ~ 2.5. We assume a ACDM cosmology with JIm = 
0.3, Qa 0.7, and Hq = 70 km s"^ Mpc"^ throughout. 
Broad-band magnitudes and colors are given in the AB 
system. 

2. DATA 

We use data from the NMBS, a moderately wide, mod- 
erately deep survey using five unique medium-bandwidth 
near-IR filters designed to sample the Balmer/4000A 
break at 1.5 < z < 3.5 with roughly twice the spec- 
tral resolution as the standard broadband JHK filters 
(van Dokkum et al. 2009). Details of the reduction 
of the medium-band filter images, object detection, and 
production of the photometric catalogs will be described 
in detail by K.E. Whitaker (in prep). Briefly, we se- 
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Fig. 1. — Data quality of the NMBS. u-8fim SEDs are shown 
for six objects at 1.7 < z < 2 that have K = 22 ± 0.1, the median 
magnitude in this redshift range. From top to bottom, the objects 
shown have Ay = [0.0,0.6,0.9, 1.5,2.1,2.5] ([HJ. The arrows point 
towards the measured MIPS fluxes and show that the objects whose 
SED fits indicate large amounts of dust are also detected at 24 fira. 



lect galaxies with K < 22.8 from two ^0.25 deg^ fields 
within the larger COSMOS and AEGIS surveys. At opti- 
cal wavelengths we use public ugriz images in both fields 
from the CFHT Legacy Survey* that were reduced by 
the CARS team (Hildebrandt et al. 2009). We include 
deep Subaru images in the B jVjr^ broad-band fil- 
ters that cover the COSMOS field (Capak et al. 2007). 
NIR images in the Spzteer-IRAC bands are provided by 
the S-COSMOS survey (Sanders et al. 2007) and over a 
strip that overlaps with ^60% of our AEGIS pointing by 
Barmby et al. (2008). Example SEDs showing the data 
quality of the NMBS are shown in Figure [1] 

We measure fluxes at 24 jxm from Spitzer-MlVS images 
provided by the S-COSMOS and FIDEL^ surveys using 
a source-fitting routine that is designed to address confu- 
sion issues in the broad-PSF MIPS images (Labbe et al. 
2006; see Wuyts et al. 2007 for an illustrative example). 
Here we are primarily concerned with whether or not an 
object is detected with 24 ^m flux, 5*24 > 20 ^Jy. We 
restrict the sample described in this Letter to the 24,800 
objects that lie within the ^0.4 deg^ overlap area covered 
by both the NMBS and MIPS images. 

We estimate photometric redshifts from the full u-8/j,m 
SEDs described above using the EAZY code (Brammer, 
van Dokkum, & Coppi 2008). We find excellent agree- 
ment between the EAZY Zphot and spectroscopic red- 
shifts measured by the zCOSMOS (Lilly et al. 2007) 

8 http://www.cfht.hawaii.edu/Science/CFHTLS/ 

9 http://irsa.ipac.caltech.edu/data/SPITZER/FIDEL/ 



and DEEP2 (Davis et al. 2003) surveys: for 632 objects 
with Zspoc < 1 from zCOSMOS, aj{l + z) = 0.016, 
and for 2316 objects with Zgpcc < 1-5 from DEEP2, 
0-2/(1 + z) = 0.018, with 3% 5f7 outliers. It is impor- 
tant to keep in mind, however, that these spectroscopic 
samples have significantly different selection functions, 
and thus will not be representative of the full X-selected 
sample used here, van Dokkum et al. (2009) targeted 
galaxies from the i^T-selected Kriek et al. (2008) sample 
to test the (NIR) medium band technique on red galax- 
ies at z > 1.7. They find crz/(l + z) = 0.01 for the four 
brightest galaxies. When 10 galaxies with S/N = 6-8 
in the NMBS filters are included the scatter increases 
to ^-^0.02, which is much smaller than has been achieved 
with broad-band filters at these redshifts. This redshift 
precision is similar to that of the COMBO- 17 survey, 
which pioneered the use of medium-band optical filters 
to determine accurate photometric redshifts and colors 
at z < 1 (Wolf et al. 2003; B04). 

We compute rest-frame U — V colors from the best-fit 
EAZY template in a similar way as done by Wolf et al. 
(2003) for COMBO-17. We find these direct template 
fluxes to be more robust than rest-frame fluxes interpo- 
lated between pairs of observed bands (i.e. Rudnick et al. 
2003; Taylor et al. 2009b) when closely-spaced medium- 
band or overlapping broad-band observed filters are used 
(Brammer, in prep). The results below remain the same 
for both methods of estimating rest-frame colors. With 
z fixed to the EAZY output, we use the SED modeling 
code, FAST (Kriek et al. 2009), to estimate the steUar 
masses, star formation rates, and dust content, parame- 
terized by the extinction in the V band {Ay), of all galax- 
ies in our sample. With FAST, we use a grid of Maras- 
ton (2005) models computed with a Kroupa (2001) IMF 
that have exponentially-declining star-formation histo- 
ries with decline rates, logr/yr = 7-10, and we allow 
dust extinction with Ay = 0-4 following the Calzetti 
et al. (2000) extinction law. Adopting an IMF without 
a turnover at low masses (e.g. Salpeter 1955) will cause 
a shift in the derived stellar masses without significantly 
changing the other properties derived from the SED fit. 

3. REST-FRAME COLOR DISTRIBUTION 

The left panel of Figure [2] shows the rest-frame U — V 
color plotted versus redshift, z, for all galaxies in our 
sample. The slope of the color-magnitude relation 
(CMR) at a given redshift was removed using the non- 
evolving OMR slope of B04, 

([/- y )' = ([/- F) -f 0.08 (My -f 20). (1) 

This correction is only valid for galaxies on the red se- 
quence and is the cause of the strong color evolution 
of the blue cloud at z < 0.5, where galaxies with low 
luminosities and correspondingly large color corrections 
enter the flux-limited sample. Our conclusions do not 
critically depend on the assumed (non-) evolution of the 
CMR slope; we flnd the same general trends described 
below assuming the evolving slope of Stott et al. (2009) 
(see Figure [3|). 

A number of observations can be made immediately 
from Figure [2l First, and a key result of this Letter 
is that a bimodal color distribution persists up to the 
highest redshifts (z ~ 2.5) where the NMBS is reason- 
ably complete. The number of red objects in the NMBS 
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Fig. 2.— Rest-frame (U - V)' color as a function of redshift for objects in both NMBS fields with K < 22.8. The color plotted is 
corrected for the slope of the (local) color-magnitude relation (Eq. [TJ. The density of points is shown as a 2D histogram scaled so that 
the maximum bin in each column is shaded black. The solid diagonal line in each panel shows the evolution of the CMR zeropoint found 
by B04. Objects detected and not detected with MIPS at 24 fira are shown in the center and right panels, respectively. Removing the 
MlPS-detected objects allows a remarkably clean separation of the red and blue galaxy sequences. The hatched regions indicate where 
our K < 22.8 sample is less than 90% complete at K'h = lO^^'^Mo. The curvature of the completeness line is due to the observed K 
magnitude selection. As discussed extensively by, e.g., Taylor et al. (2009a), ii'-selected samples are more complete for blue galaxies than 
for red galaxies at a fixed mass. 



drops off sharply at z 2, which is hkely due to in- 
completeness for all but the highest stellar masses and 
a consequence of the K-hand selection. Proper analysis 
of the number density evolution of the red and blue se- 
quences will require more careful selection based on more 
physical properties, such as stellar mass, which we will 
present in a forthcoming paper (Brammer, in prep). A 
second result clear from Figure [2] is that the color of the 
red peak evolves slowly with redshift and is a natural 
extension of the red peak found at < z < 1 by B04. 
The color evolution of the blue peak is rather similar. 
Finally, the MlPS-detected objects have rest-frame col- 
ors that fall between the red and blue sequences (see also 
Figure [3|), which has been found previously by Cowie & 
Barger (2008) and Salim et al. (2009) at z < 1.5 but 
that we extend here to significantly higher redshifts (see 
also W09). The contrast between the red and blue se- 
quences is greatly enhanced when the MlPS-detected ob- 
jects are removed. Note that here we are considering as 
the red sequence the red peak of CMR-corrected colors 
as a function of redshift. We show that this interpreta- 
tion is consistent with the more traditional definition of 
the red sequence in the color-mass diagram in 53 

4. DUST IN THE GREEN VALLEY 

The U — V colors of the MlPS-detected galaxies sug- 
gest that the green valley is largely populated by dusty 
star-forming galaxies. Figure [3] shows Ay from the SED 
fit as a function of U — V color, which is now corrected 
for both the slope of the CMR and for the B04 evolution 
of the CMR zeropoint with redshift (see Figure [2]). Ob- 
jects on the red sequence are clearly visible as the clump 
at corrected U — V ^ 2 with low Ay values and very low 
sSFR. Some galaxies in the red clump have Ay — 0.5 or 
larger, which is perhaps larger than might be expected 
for truly "red and dead" galaxies, but the la error on 
Ay estimated by FAST includes zero for most of these 
galaxies. The median Ay (thick black curve in Figure[3]) 
peaks in the green valley, demonstrating that it is indeed 
populated at redder colors by increasingly dusty objects 
(see also W07; W09). A green valley with colors deter- 
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Fig. 3. — Ay vs. rest-frame U — V color [corrected for the 
CMR slope and its evolution with redshift (B04; see Figure [2j] for 
galaxies at 1.0 < z < 2.5. The points are color-coded by the specific 
SFR (sSFR) determined from the SED fit. A grid of AAy = 0.1 
was used for the fit; here the points are shown with a small random 
offset for clarity. The thick black line and gray shaded region show 
the median and interquartile range of Ay as a function of color. 
The dotted line is the median trend for colors corrected using the 
Stott et al. (2009) CMR slope evolution. The reddening vector is 
shown for the Calzetti et al. (2000) law with AAv = 1. The color 
distributions of all MlPS-detected objects and of objects above and 
below Ay = 1 are shown in the bottom panel. The vertical dashed 
lino shows the red-sequence color selection used by B04. 

mined primarily by aging stellar populations would be 
horizontal on this figure. 

The bottom panel of Figure [3] (see also Figure [T]) shows 
that the color distribution of objects with MIPS detec- 
tions (negative red histogram) is similar to the distri- 
bution of objects with Ay > 1 (positive line-filled his- 
togram). The relative contributions of star- formation 
and AGN to the 24 /im flux are difficult to assess here. 
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Fig. 4. — The color-mass relation at 1.8 < z < 2.2. The U ~ V color in the right panel is corrected for reddening with the factor from 
Eq. [21 Objects detected at 24 /^m are shown as filled gray circles. The hatched region indicates the approximate 90% completeness limit 
for red-sequence objects at z = 2. Due to the reddening correction, the completeness limit for blue galaxies is more difficult to estimate. 
The colorbar in the right panel shows the median sSFR as a function of dust-corrected color, with the same coding as Figure [3] 



particularly because the adopted MIPS flux limit corre- 
sponds to a large range of MIR luminosity, sampled at a 
broad range of rest-frame wavelengths, over the redshift 
range considered. Only ^1% of the galaxies in Figure 
[3] have X-ray counterparts, which implies that powerful 
AGN are probably not the dominant source of MIPS flux 
in the green valley. 

5. THE DUST-CORRECTED RED SEQUENCE AT 

Z = 2 

Figure [3] demonstrates that the overlapping dusty- 
starburst and passive galaxy populations at the reddest 
U — V colors have SEDs that can be easily distinguished 
with the NMBS photometry. We explore this separation 
further in Figured! which shows the color-mass diagram 
at 1.8 < z < 2.2 where the the right panel shows the 
rest-frame U — V color corrected for dust reddening 
with the correction factor, 

Auv=0A7-Av, (2) 

computed for the Calzetti et al. (2000) extinction law. 
After the reddening correction the galaxies collapse into 
strikingly well-deflned red and blue sequences, now in the 
traditional sense of the terms. We remark, though, that 
the color-mass diagram shown is somewhat more model- 
dependent than a color-magnitude diagram. Figure |4] 
extends to signiflcantly higher redshift and with good 
statistics the results found by Wyder et al. (2007) {z ~ 
0) and Cowie & Barger (2008) {z ^ 1.5), who show that 
the blue and red sequences are more easily separated in 
the CMD if one considers de-reddened colors. Similarly, 
many red objects at z ^ have intrinsic colors that put 
them on the blue sequence when the extinction effects 
of inclined disk galaxies are taken into account (Martin 
et al. 2007; Baihn & Harris 2008; Mailer et al. 2009). 

There is some question as to how well Ay can be es- 
timated from the SED fit, since there are complicated 
degeneracies between the parameters of the fit (Ay, age, 
SFH, Z). Figure [T] demonstrates that there are real 



differences in the SED shapes for galaxies with simi- 
lar [/ — V colors but different values of derived Ay. 
The combination of the medium-band NIR filters and 
the IRAC photometry at Arost > 1/im traces the detailed 
shape of the SED near the Balmer/4000A break and de- 
termines the SED slope redward of the break (i.e. the 
V — J color), which allows the separation of dusty and 
passive SED types with similar U ~ V colors (Labbe 
et al. 2005; W07; W09). 

6. DISCUSSION 

We have shown that either removing MlPS-detected 
galaxies or applying a reddening correction to rest-frame 
U — V colors allows us to cleanly divide the pop- 
ulation of galaxies with red rest-frame colors among 
dusty-star burst and "red and dead" galaxies. The lat- 
ter, which have intrinsically red colors characteristic of 
evolved stellar populations, form what we call the "dead 
sequence", which is in place at all < z < 2.5. This 
term was used previously by Romeo et al. (2008), who 
consider the sub-sample of red-sequence galaxies with 
log(sSFR-yr) < -11; Figure S] shows that the cloud of 
galaxies with red dust-corrected colors do, in fact, have 
very low sSFR. We note, however, that there could be 
galaxies on what we consider the "dead sequence" with 
somewhat higher sSFR than in Romeo et al. (2008). 

Obtaining these results at z > 1.5 is made possible 
by the NEWFIRM Medium-band Survey, which pro- 
vides high-quality Zphot and well-sampled SEDs of red 
galaxies at redshifts beyond the current limits of large 
spectroscopic surveys (van Dokkum et al. 2009): of the 
'^SOOO NMBS objects with measured spectroscopic red- 
shifts, only 9 have U — V > 1.4 at z > 1.2. Conversely, 
there are 2826 objects that satisfy these criteria in the 
full photometric sample. 

The results presented here are consistent with those 
of W09, who show that quiescent galaxies can be sepa- 
rated from dusty star-forming galaxies in the U — V 
vs. V — J color-color diagram out to z — 2. The 
reddening-corrected dead sequence galaxies discussed 
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here fall nicely within the UVJ quiescent selection cri- 
teria, though the reddening correction is critical for de- 
tecting the bimodal color distribution in a color-mass 
diagram at z > 1.5, explaining why it is not seen by 
W09. 

We find the color evolution of dead sequence galaxies 
to be consistent with that found by B04, who show that 
it can be roughly explained by simple stellar populations 
formed at z > 2 that age passively to the present day. We 
show, however, that additional care must be taken when 
defining red sequence galaxy samples as done by, e.g., 
B04. The bottom panel of Figure [3] shows that a red 
galaxy selection based on a single rest-frame color will 
contain a significant fraction of dusty-starburst galax- 
ies (~20% with Ay > 1) that would otherwise be much 
bluer (see also Cimatti et al. 2002; Brammer & van 
Dokkum 2007; W07; Gallazzi et al. 2009; W09; Wolf 
et al. 2009). The reddening correction dramatically re- 
duces the number of objects in the green valley, so such a 
correction would be necessary when using those objects 
to estimate the rate at which galaxies move from the 
blue to red sequences (Martin et al. 2007). For exam- 
ple, a sparsely-populated green valley at most redshifts 
could constrain the rate at which star-forming galaxies 
are quenched, which in turn may help constrain the pro- 
cesses that regulate star formation in theoretical models. 

As might be expected, the large majority of objects 
detected at 24 ^m have reddening-corrected colors that 
place them on the blue-sequence (Reddy et al. 2006, see 
also Figure 3). A number of MlPS-detected galaxies, 
however, remain on the red sequence even after apply- 
ing the reddening correction: roughly one in seven MIPS 
objects have dust-corrected colors within 0.5 mag of the 
red sequence ridgeline drawn in Figure [21 The dead se- 
quence galaxies detected with MIPS have de-reddened 
colors that are slightly bluer than those without MIPS 
detections (Figured]). This could perhaps be due to low 
levels of recent star formation or to AGN contamination 
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of the UV-optical light, though this tentative result re- 
quires more detailed investigation because, for example, 
it is not clear that the dust correction or even the extinc- 
tion law itself are appropriate for these galaxies. 

Figures [2] and [4] demonstrate that the dead sequence 
persists to z = 2-2.5, which is consistent with the dis- 
covery of a small (spectroscopic) sample of red sequence 
objects with low sSFR at z ~ 2.3 described by Kriek 
et al. (2008). Such studies of red sequence galaxies found 
among the field are complementary to others that target 
the higher density environments of clusters at z > 1 (e.g., 
Mei et al. 2009) and proto-clusters at z ~ 2 (Zirm et al. 
2008). Mechanisms that have been proposed to suppress 
star formation in massive galaxies, including the much- 
discussed feedback from AGN and gravitational heating 
from cosmological accretion (Naab et al. 2007; Dekel 
& Birnboim 2008), must be able to explain the forma- 
tion of the dead sequence in a variety of environments 
as early as z ^ 2.5. Finally, we note that quantitative 
constraints on the assembly and evolution of passively- 
evolving galaxies require studies of the evolution of their 
mass function, possibly combined with structural infor- 
mation (see, e.g., B04; Faber et al. 2007; Cimatti et al. 
2008; van Dokkum et al. 2008, and many other stud- 
ies). The NMBS is ideally suited to help disentangle 
these processes by providing accurate photometric red- 
shifts and colors over a large area in the redshift regime 
of the most rapid build-up of massive galaxies. 
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